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Abstract A major contribution to the quest of constructing quantum dynamics of non-Abelian 
fieds is due to V.N. Gribov. Perturbative approach to the colour confinement, both in gluo- 
dynamics and the real world, was long considered heretic but is gaining ground. We discuss 
Gribov's approach to the confinement problem, centered around the role played by light quarks 
— the supercritical light quark confinement scenario. We also review some recent developments 
that are motivated, directly or indirectly, by his ideas. 
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1 INTRODUCTION 

One of the most challenging problems for modern theoretical physics is under- 
standing the confinement of colour — the selection rule suggested by quantum 
chromodynamics (QCD) to explain the fact that quarks (and gluons) appear 
in the physical spectrum only imprisoned inside composite states — colourless 
mesons and baryons. 

An amazing success of the relativistic theory of electron and photon fields - 
quantum electrodynamics (QED) — has produced a long-lasting negative impact: 
it taught the generations of physicists that came into the business in/after the 
70's to "not to worry" . Indeed, today one takes a lot of things for granted. 

One rarely questions whether the alternative roads — secondary quantization, 
functional integral and the Feynman diagram approach — really lead to the same 
quantum theory of interacting fields. 

One feels ashamed to doubt an elegant and powerful, but potentially deceiving, 
technology of translating the dynamics of quantum fields into that of statistical 
systems (imaginary time = Euclidean rotation trick, resulting in the coupling <-> 
temperature analogy). 

One takes the original concept of the "Dirac sea" — the picture of the fermionic 
content of the vacuum — as an anachronistic model, a sort of Maxwell's mechan- 
ical ether, Fd-rather- do- without. 

One was taught to consider the problems with field-theoretical description of 
point-like objects and their interactions at very small distances (ultraviolet di- 
vergences) as purely technical: renormalize it and forget it. 

So far so good for QED, where the physical objects — electrons and photons 
- are direct images of the fundamental fields that one put into the local La- 
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grangian of the theory. In other words, interacting fields closely resemble their 

bare counterparts. In these circumstances, the ground state of the theory — the 
physical vacuum — turns out to be, in a certain sense, trivial. The role played 
by the vacuum fluctuations of electromagnetic fields and by the "Dirac sea" of 
negative-energy electron states, is as important as it is straightforward. Interac- 
tion with the vacuum makes the effective interaction strength a e . m . (as well as 
the effective electron mass) run with the distance. At the same time, it does not 
affect the nature of the interacting fields. 

Not so clear in QCD. Here the physical hadronic states are not in one-to-one 
correspondence with the fundamental quarks and gluons: interaction with the 
QCD vacuum changes the bare fields beyond recognition. In such an environment, 
our preconceived ideas about QFT dynamics may turn out to be a handicap, a 
serious hidden obstacle on the way to attacking the confinement problem. 

Vladimir Naumovich Gribov (1930-1997) belonged to a generation of physicists 
(now almost extinct) that thought about the Quantum Field Theory (QFT), that 
was used to questioning its foundations, concepts and means. An invitation to 
the "Gribov conception of QCD" is an invitation to unlearn. Learning to unlearn 
isn't easy. This possibly explains why the programme that Gribov formulated and 
was pursuing in the 90's of explaining the confinement of colour as "supercritical 
binding" of light quarks has yet to receive the attention it deserves from the 
physics community at large. 

It goes without saying that Gribov could be wrong in his vision about the 
nature of the QCD confinement. However, it does not seem to us a good policy 
to be simply indifferent to what one of the creators of the modern theoretical 
physics had to say on the subject, to the ideas and tools he has developed during 
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the last half of his professional life. 

This review does not attempt to cover Gribov's impact on theoretical physics. 
Suffices it to mention that his name is attached to many a key notion of the the- 
ory arsenal: Gribov-Froissart projection and the Gribov vacuum pole (Pomeron), 
Gribov factorization, Reggeon Calculus and Reggeon Field Theory (RFT), Gri- 
bov diffusion, the AGK (Abramovsky-Gribov-Kancheli) cutting rules, the Gribov 
bremsstrahlung theorem, Gribov-Glauber theory of relativistic multiple scatter- 
ing, Gribov-Lipatov evolution equations, Gribov copies and the horizon, etc. 

Let us remark that Gribov's impact on modern physics is deeper than it is 
generally known to be. A couple of examples will illustrate the point. 

Working on the problem of the so-called strong-coupling regime of interacting 
Pomerons, V. Gribov and A. Migdal developed an ingenious technique for ana- 
lysing dynamical systems with long-range fluctuations. Such fluctuations being 
typical for condensed matter physics near the critical temperature, this trig- 
gered an important breakthrough in solid state physics. The subsequent works 
of A. Migdal and A. Polyakov, and a contemporary more general treatment sug- 
gested by L. Kadanoff and K. Wilson, have established the scaling solution of the 
problem of the second order phase transitions. 

In 1969, in one of his jewels Interaction of photons and electrons with nuclei at 

high energies Gribov established and described the space-time picture of particle 

interactions at high energies. This work found a way through the iron curtain. 1 

Its key elements were incorporated into the famous Feynman book which laid the 

foundation of the parton model. The Feynman-Gribov parton model, that is. 

The last 20 years of his life V. Gribov devoted to non-Abelian quantum gauge 
x It was reported, prior to publication, by J. Bjorken at Caltech seminar. 



6 Dokshitzer & Kharzeev 

theories, and to the QCD confinement problem in particular. His very first QCD 

study of 1976-77 produced a brilliant physical explanation of asymptotic freedom 
introduced into the physics of strong interactions by D. Gross, F. Wilczek and 
H. Politzer. Gribov's approach was based on an early observation of the anti- 
screening phenomenon made by I. Khriplovich in a pre-historic 1969, and revealed 
the inconsistency of the standard field-theoretical treatment of gluon fields (Gri- 
bov copies). In the late 80's Gribov suggested the quark confinement scenario 
based on the so-called supercritical binding of light fermions by a quasi-Coulomb 
interaction. 

His last works remain to be discovered, understood and developed. 

For V.N. Gribov, the problem of confinement as the problem of understanding 
the dynamics of vacuum fluctuations, and of the structure of hadrons as physical 
states of the theory, was always inseparable from the problem of understanding 
the physics of high energy hadron scattering (the Pomeron picture) . 

Therefore, we start by an overview of the Gribov's development of high-energy 
physics. The two subsequent Sections are devoted to his contribution to the 
quantum dynamics of non-Abelian Yang-Mills gauge fields in general, and to the 
basics of the Gribov light-quark confinement picture. The last Section briefly de- 
scribes the byproducts of the QCD programme: an unpublished study of QED at 
short distances (a potential resolution of the notorious "Landau ghost" problem) 
and his picture of a composite (super-critically bound) Higgs boson. 

2 HADRON INTERACTIONS AT HIGH ENERGIES 

In the late fifties, when Vladimir Gribov, then a young researcher at the Ioffe 
Physico- Technical Institute (Leningrad, USSR) became interested in the physics 
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of strong hadron interactions, not much was understood about processes at high 
energies. The only theoretical result, derived from the first principles, was the 
Pomeranchuk theorem — an asymptotic equality of particle and antiparticle total 
scattering cross sections. 

2.1 Asymptotic Behaviour s a w 

Gribov's 1960 paper Asymptotic behaviour of the scattering amplitude at high 
energies (^Q) was a breakthrough. Using the so-called double-dispersion represen- 
tation for the scattering amplitude, suggested by S. Mandelstam back in 1958 ©, 
Gribov proved an inconsistency of the then popular black disk model of hadron- 
hadron scattering. This analysis may be considered as the first building block put 
into the edifice of the modern theory of hadron interactions. It has demonstrated 
the combined power of the general principles of relativistic quantum theory - 
unitarity (conservation of probability), analyticity (causality) and the relativistic 
nature (crossing symmetry) — as applied to high-energy scattering phenomena. 

By studying the analytic properties in the cross channels, he showed that the 
imaginary part of the scattering amplitude in the form 



that constituted the black disk model of diffraction in the physical region of 
s-channel scattering, contradicts the unitarity relation for partial waves in the 
crossing t-channel. To solve the puzzle, Gribov suggested the behaviour of the 
amplitude (for large s and finite t) in the general form 



where Bt is a slow (non-exponential) function of In s (decreasing fast with t) and 



A 1 (s,t) = sf(t) 



(1) 




(2) 
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q(0) = 1 ensures the approximate constancy of the total cross section with energy. 

In this first paper Gribov only analysed the constant exponent, q(t) = 1, 
having remarked on the possibility q(t) / const as "extremely unlikely" . Indeed, 
considering the t-dependence of the scattering amplitude, this would correspond 
to a strange picture of the radius of a hadron infinitely increasing with energy. 2 
(In this particular case he proved that the cross section has to decrease at high 
energies, Bt(lns) < 1/lns, to be consistent with the t-channel unitarity.) 

Soon after that Gribov became aware of the finding by T. Regge © that in 
non-relativistic quantum mechanics, in the unphysical region \t\ 3> s (momentum 
transfer much larger than the energy, corresponding to large imaginary scattering 
angles | cos@| — * oo), the scattering amplitude has a form 



where £(s) is the position of the pole of the partial wave f# in the complex plane 
of the orbital momentum i? 

T. Regge found that the poles of the amplitude in the complex £-plane were 
intimately related with bound states/resonances. It is this aspect of the Regge 
behaviour that initially attracted the most attention: 

"5". Mandelstam has suggested and emphasized repeatedly since 1960 
that the Regge behaviour would permit a simple description of dynam- 
ical states (private discussions). Similar remarks have been made by 

R. Blankenbecker and M.L. Goldberger and by K. Wilson." (@J 
2 Reportedly, it was L.D. Landau who "forced" Gribov to publish Eq. J2J in its general form, 

with q(t) 7^ const. "You are too young to judge" the blessing went, according to legend. 

3 Gribov apparently learned about the Regge result from a paper by G. Chew and S. Frautschi 

of 1960 which contained a footnote describing the main Regge findings. 



A(s,t) oc t^ s) 



(3) 
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The structure of the Regge amplitude Eq. @ motivated Gribov to return to 

the consideration of the case of i-dependent exponent in his general high-energy 

ansatz Eq. (J2J) that was dictated by t-channel unitarity. His letter to ZhETF 

Possible asymptotic behaviour of elastic scattering (5 ) became the first application 

of Regge ideas to the high-energy asymptotic behaviour of scattering amplitudes. 

By then M. Froissart had already proved his famous theorem that limits the 

asymptotic behaviour of the total cross sections ©, 

<r tot oc s- 1 0)1 < C ln 2 s. (4) 

Thus, having accepted £(0) = 1 for the rightmost pole in the £-plane — the vac- 
uum pole — as the condition "that the strongest possible interaction is realized" , 
Gribov formulated and discussed "the main properties of such an asymptotic scat- 
tering behaviour . . . which in spite of having a few unusual features is theoretically 
feasible and does not contradict the experimental data": 

• the total interaction cross section is constant at high energies, 

• the elastic cross section tends to zero as 1/lns, 

• the scattering amplitude is essentially imaginary, 

• the significant region of momentum transfer in elastic scattering shrinks 
with energy increasing, y/—i oc (Ins) -1 / 2 . 

He also analysed the s-channel partial waves to show that for small impact pa- 
rameters p < R their amplitudes fall as 1/lns, while the interaction radius R 
increases with energy as p oc vTn s. He concluded: 

"this behaviour means that the particles become grey with respect to 
high-energy interaction, but increase in size, so that the total cross 
section remains constant." 
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Thus, shrinkage of the diffractive peak was predicted, which was then experimen- 
tally verified at particle accelerator experiments in Russia (IHEP, Serpukhov), 
Switzerland (CERN) and the US (FNAL, Chicago), as were the general relations 
between the cross sections of different processes, that followed from the Gribov 
factorization theorem (JJJ). 

These were the key features of what has become known, quite ironically, as the 
"Regge theory" of strong interactions at high energies. 

On the opposite side of the curtain, the basic properties of the Regge pole pic- 
ture of forward/backward scattering were formulated half a year later by G. Chew 
and S. Frautschi (|5|). In particular, they suggested "the possibility that the re- 
cently discovered p meson is associated with a Regge pole whose internal quantum 
numbers are those of an I = 1 two-pion configuration" and conjectured the uni- 
versal high-energy behaviour of backward 7r + 7r°, K + K° and pn scattering due 
to p-reggeon exchange. Chew and Frautschi also stressed that the hypothetical 
Regge pole with a(0) = 1 responsible for forward scattering possesses quantum 
numbers of the vacuum. 

Dominance of the Gribov vacuum pole automatically validates the Pomer- 
anchuk theorem. The name "Pomeron" for the vacuum pole was coined by Mur- 
ray Gell-Mann, who referred to Geoffrey Chew as an inventor. 

2.2 Complex Angular Momenta in Relativistic Theory 

In non-relativistic quantum mechanics the interaction Hamiltonian allows for 
scattering partial waves to be considered as analytic functions of complex an- 
gular momentum I. Gribov's paper Partial waves with complex orbital angular 
momenta and the asymptotic behaviour of the scattering amplitude showed that 
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the partial waves with complex angular momenta can be introduced in a relativis- 

tic theory as well. Here it is the unitarity in the crossing channel that replaces the 

Hamiltonian and leads to analyticity of the partial waves in I, The corresponding 

construction is known as the "Gribov-Froissart projection" Q. 

Few months later Gribov demonstrated that the simplest (two-particle) t- 
channel unitarity condition indeed generates moving poles in the complex ^-plane. 
This was the proof of the Regge hypothesis in relativistic theory (|1U[). 

The "Regge trajectories" a{t) combine hadrons into families: Sh = a(rnV), 
where Sh and mh are the spin and the mass of a hadron (hadronic resonance) with 
given quantum numbers (baryon number, isotopic spin, strangeness, etc.) (JHJ). 
Moreover, at negative values of i, that is in the physical region of the s-channel, 
the very same function a(t) determines the scattering amplitude, according to 
Eq. (J2J). It looks as if high-energy scattering was due to i-channel exchange of a 
"particle" with spin a(t) that varies with momentum transfer t — the "reggeon". 

Thus, the high-energy behaviour of the scattering process a + b — ► c + dis 
linked with the spectrum of excitations (resonances) of low-energy scattering in 
the dual channel, a + c — > b + d. This intriguing relation triggered many new ideas 
(bootstrap, the concept of duality). Backed by the mysterious linearity of Regge 
trajectories relating spins and squared masses of observed hadrons, the duality 
ideas, via the famous Veneziano amplitude, gave rise to the concept of hadronic 
strings and to development of string theories in general. 

2.3 Interacting Pomerons 

A lot of theoretical effort was invested into understanding of the approximately 
constant behaviour of total cross sections at high energies. 
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To construct a full theory that would include the Pomeron trajectory with 

the maximal "intercept" ap (0) = 1 that respects the Froissart bound, and would 

be consistent with unitarity and analyticity proved to be very difficult. This is 

because multi- Pomeron exchanges become essential, which generate branch points 

in the complex £-plane, which singularities were first discovered by Mandelstam 

in his seminal paper (|11[). Moreover, the study of particle production processes 

with large rapidity gaps led Gribov, Pomeranchuk and Ter-Martirosyan to the 

concept of interacting reggeons. By the end of the 60-s V. Gribov had proposed 

the general theory known as Gribov Reggeon Calculus. He formulated the rules 

for constructing the field theory of interacting Pomerons — the Reggeon Field 

Theory (RFT) — and developed the corresponding diagram technique. Gribov 

RFT reduces the problem of high energy scattering to a non-relativistic QFT of 

interacting particles in 2+1 dimensions. 

2.4 Gribov Partons and Feynman Partons 

One of Gribov's most important contributions to high energy hadron physics 
was the understanding of the space-time evolution of high energy hadron-hadron 
and lepton-hadron processes, in particular the nature of the reggeon exchange 
from the s-channel point of view. In 1973, in his lecture at the LNPI Winter 
School (|12|). Gribov outlined the general phenomena and typical features that 
were characteristic for high energy processes in any QFT. 

To understand the structure of hard (deep inelastic) photon-hadron interac- 
tions Feynman suggested the idea of partons — point-like constituents of hadronic 
matter. Feynman defined partons in the infinite momentum frame to suppress 
vacuum fluctuations whose presence would have undermined the notion of the 
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parton wave function of a hadron (|13f) . 

The power of Gribov's approach lied in applying the universal picture of fluc- 
tuating hadrons to both soft and hard interactions. Gribov's partons are con- 
stituents of hadron matter, components of long-living fluctuations of the hadron 
projectile, which are responsible for soft hadron-hadron interactions: total cross 
sections, diffraction, multiparticle production, etc. 

Gribov's earlier work Interaction of j-quanta and electrons with nuclei at high 
energies Q14|) had been a precursor to the famous Feynman paper. Gribov de- 
scribed the photon interaction in the rest frame of the target nucleus. An incident 
real photon (or a virtual photon in the deep inelastic scattering case) fluctuates 
into hadrons before the target, at the longitudinal distance L increasing with 
energy. 4 Therefore, at sufficiently large energy, when the fluctuation distance ex- 
ceeds the size of the target, the photon no longer behaves as a point-like weekly 
interacting particle. Its interaction resembles that of a hadron and becomes 
"black" , corresponding to complete absorption on a large nucleus. 

Being formally equivalent to Feynman's treatment, Gribov's approach is better 

suited for the analysis of deep inelastic phenomena at very small Bjorken x, where 

the interaction becomes actually strong, and the perturbative QCD treatment is 

bound to fail. Gribov diffusion in the impact parameter space giving rise to 

energy increase of the interaction radius and to the reggeon exchange amplitude, 

coexisting fluctuations as a source of branch cuts, duality between hadrons and 

partons, a common basis for hard and soft elastic, diffractive and inelastic process 

- these are some of the key features of high energy phenomena in QFTs, which 
4 For the e~p deep inelastic scattering B.L. Ioffe has shown that the assumption of Bjorken 

scaling implies L ~ l/2xmjv, with x the Bjorken variable and mjv the nucleon mass GEJ. 
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are still too hard a nut for QCD to crack. 
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2.5 Gribov Reggeon Field Theory 

The two best known applications of the Gribov RFT are 

• general quantitative relation between the shadowing phenomenon in hadron- 
hadron scattering, the cross section of diffractive processes and inelastic 
multi-particle production, known as "Abramovsky-Gribov-Kanchelli cut- 
ting rules" (AGK) JTSJ), and 

• the essential revision by Gribov of the Glauber theory of nuclear shadowing 
in hadron- nucleus interactions (17). 

In 1968 V.N. Gribov and A. A. Migdal demonstrated that the scaling behaviour 
of the Green functions emerged in QFT in the strong coupling regime (|18|). As 
we have mentioned in the Introduction, their technique helped to build the quan- 
titative theory of second order phase transitions and to analyse critical indices 
characterising the long-range fluctuations near the critical point. 

The problem of high energy behaviour of soft interactions remained unsolved, 
although some viable options were suggested. In particular, in Properties of 
Pomeranchuk Poles, diffraction scattering and asymptotic equality of total cross 
sections (19 ) Gribov has shown that a possible consistent solution of the RFT in 
the so-called weak coupling regime calls for the formal asymptotic equality of all 
total cross sections of strongly interacting particles. 

Gribov's last work in this subject was devoted to the intermediate energy range 
and dealt with interacting hadron fluctuations (20). 

The study of the strong coupling regime of interacting reggeons (pioneered by 
A.B. Kaidalov and K.A. Ter-Martirosyan) led to the introduction of the bare 
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Pomeron with ap(0) > 1. The RFT based on i-channel unitarity should enforce 

the s-channel unitarity as well. The combination of increasing interaction radius 

and the amplitudes in the impact parameter space which did not fall as 1/lns (as 

in the one-Pomeron picture) led to logarithmically increasing asymptotic cross 

sections, resembling the Froissart regime (and respecting the Froissart bound 

Eq. ©). The popularity of the notion of the bare Pomeron with ap(0) > 1 is 

based on experiment (increasing Ctot)- Psychologically, it is also supported by 

the perturbative QCD finding that the (small) scattering cross section of two 

small-transverse-size objects should increase with energy in a power-like fashion 

in the restricted energy range — the famous "hard BFKL Pomeron" (|21j) . 

2.6 Reggeization and Pomeron Singularity in Gauge Theories 

In the mid-sixties Gribov initiated the study of double logarithmic asymptotics of 
various processes in QED, making use of the powerful technique he had developed 
for the analysis of the asymptotic behaviour of Feynman diagrams (22). 

In particular, in 1975 Gribov, Lipatov and Frolov studied the high energy 
behaviour of QED processes from the point of view of "Regge theory". High 
energy scattering amplitudes with exchange of an electron in the i-channel ac- 
quire, in higher orders in QED coupling, a characteristic behaviour A cx s J W with 
j( m e) = 1/2- This means that electron becomes a part of the Regge trajectory: 
reggeizes. 5 For the vacuum channel, however, they found (23) that the rightmost 
singularity in the complex j-plane is not a moving pole (as it is for electron) 
but, instead, a fixed branch point singularity positioned to the right from unity, 

j(0) = 1 + co? > 1. This was a precursor of a similar result found later by Fadin, 
5 And so do quarks and gluons in QCD; Fadin, Frankfurt, Lipatov, Sherman (1976). 
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Lipatov and Kuraev (121[) in non-Abelian theories, and QCD in particular. The 

problem of apparent anti-Froissart behaviour of the perturbative "hard Pomeron" 
in QCD still awaits resolution. 

With the advent of QCD as a microscopic theory of hadrons and their in- 
teractions, the focus of theoretical studies has temporarily shifted away from 
Gribov-Regge problematics to "hard" small-distance phenomena. 

3 NON-ABELIAN GAUGE THEORIES 

V.N. Gribov became interested in non-Abelian field theories in 1976. His very 
first study, as a QCD apprentice, produced amazing results. In February 1977 in 
the proceedings of the 12 th Leningrad Winter School he published two lectures 
which were to change forever the non-Abelian landscape (|24U25|) . 

3.1 Anatomy of Asymptotic Freedom 

In the first lecture Instability of non-Abelian gauge fields and impossibility of the 
choice of the Coulomb gauge (|24|) Gribov gave an elegant physically transparent 
explanation of the asymptotic freedom introduced by D. Gross & F. Wilczek and 
H.D. Politzer in 1973 

The anti- screening phenomenon had been first observed back in 1969 (|27f) 
for the non-Abelian SU(2) theory in the ghost-free Coulomb gauge within the 
Hamiltonian approach. In the Hamiltonian language, there are (or rather seem to 
be) N 2 —l massless B quanta (transverse gluons, _L) and, as in QED, an additional 
Coulomb field (0). It is important to stress that the latter is not a physical 
quantum degree of freedom but a means for describing classical instantaneous 
interaction between colour charges. Unlike QED, the non-Abelian Coulomb field 
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has a colour charge of its own. Therefore, traversing the space between two 

external charge, it may virtually decay into two transverse fields, 

-> _L + _L -» , (5a) 

or into a qq pair 

-» q + q -» 0, (5b) 

in the same manner as the QED Coulomb field fluctuates in the vacuum into 
an e + e~ pair. Both these effects lead to screening of the colour charge of the 
external sources, in a perfect accord with one's physical intuition. 

This (anti- asymptotically- free) behaviour of the running coupling was first 
found by Landau, Abrikosov and Khalatnikov for QED (28). It turned out to 
be common for all then-known renormalizable field theories: with scalar (A</> 4 ), 
Yukawa, four-fermion interactions (J29J) as well as for pedagogically valuable Lee 
model (j3U|) . 

This observation had dramatic consequences: the physical interaction (renor- 
malized coupling) was predicted to vanish in the limit of a point-like bare in- 
teraction, Auv 00 • in the late 1950s the problem was known as "Moscow 
Zero". The depth of the crisis can be measured by the Dyson prophecy 1)31(1 
that the correct "meson" theory - the theory of strong interactions - 11 will not 
be found in the next hundred years" and by the Landau conclusion (32 ) that "the 
Hamiltonian method for strong interactions is dead and must be buried, although 
of course with deserved honour? 

Universality of the screening phenomenon was readily understood as a conse- 
quence of unitarity in the cross-channel. Indeed, the discontinuity of the loop 
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diagram at t > 0, 
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oc a 



describes, according to Cutkosky rules, production cross section of two physical 
(sic!) particles in the t-channel. Positivity of the imaginary part at t > then 
dictates the screening sign of the virtual loop in the s-channel (s > 0, t < 0). 

Thus, the fact that there are "physical" fields in the intermediate state in 
Eq. (JjJJ) — quark-antiquark or two transverse gluons — fixes the sign of the 
virtual correction to correspond to screening, via the unitarity relation in the 
cross-channel. Technically speaking, these virtual decay processes contribute to 
the QCD /3-function as 



daJ 1 (R)\ 



phys 



dlnR 



1 »r 2 

oc - Is H nt , 

3 3 ; 



(6a) 



that is, make the effective coupling decrease at large distances R between the 
external charges, as in QED. 

Where then the anti- screening comes from? It originates from another, specif- 
ically non-Abelian, effect namely, interaction of the Coulomb field with the field 
of "zero-fluctuations" of transverse gluons in the vacuum, 



+ 











+ .... 



In a course of such multiple rescattering, the Coulomb "quantum" preserves its 
identity as an instantaneous interaction mediator, and therefore is not affected 
by the unitarity constraints. For n = 1 the contribution vanishes upon averag- 
ing. Statistical average over the transverse vacuum fields in the second order of 



perturbation theory (n = 2) results in an additional contribution to the Coulomb 
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interaction energy which, translated into the running coupling language, gives 

{ w « 

According to Gribov, an anti-intuitive minus sign in Eq. (|6b|) has its own simple 
explanation. It is of the same origin as the minus sign in the shift of the energy 
of the ground state of a quantum-mechanical system under the second order in 
perturbation: 

„ = = „. 

The role of perturbation 8V is played by the vacuum field of transverse gluons. 

Taken together, the two contributions Eq. © combine into the standard QCD 
/^-function: 

d a s (Q 2 ) , fa s (Q 2 )\ 2 , „ . 3 , , _ 11JV C 2N f 



d a s (Q 2 )_ ( a s (Q 2 ) Y ( 3] 



3.2 Infrared Instability 

The three-dimensional transversality condition 

8 A ■ 

(V-A) = ^ = 0, i = l,2,3 (7) 

is usually imposed on the field potential to describe massless vector particles (the 
Coulomb gauge). Being antisymmetric, the field strength tensor F£ v does not 
contain time derivative of the zero-component of the potential Aq. In the Hamil- 
tonian language, this puts Aq in a position of a cyclic variable which does not 
constitute a physical degree of freedom. It can be eliminated from field dynam- 
ics H24U33|) . contributing to the Hamiltonian that describes transverse gluons, 



drx 



j-\ E a ± J ■ E\J~2 + B a ± ■ B a ± 
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an additional term responsible for Coulomb interaction between "charges", 



tfcoul = 3 / d 6 x d d y J-2p a (x) J K ab (x, y; A ± ) p\y) J~2 



(8b) 



where p a is given by the sum of the colour charge density of external sources 
(e.g., static quarks) and that induced by the transverse gluons themselves, 

P = Pq + ig 8 [A_l,E_l]. (8c) 

Introducing the covariant derivative D, 

D[A]C = VC + ig s [A,C], 



the Coulomb energy "propagator" K in Eq. (|8b[) reads 



K ab {x,y;A) 



ab 



(9) 



_D[A] • V D[A] • V 

We see that the propagation of a Coulomb field Aq in the "external" transverse 
gluon field A^ is governed by the operator that resembles the propagator of the 
Faddeev-Popov ghost, 



(D • V) Aq = V 2 A + ig s [A ± , VA ] . 



The factor J in Eq. (jSJ) is the determinant of this operator, 



(10) 



J = J[A] = -det(D[A] • V). 

Taking the expectation value of K(x, y; A) over the vacuum fields A^ produces 
an instantaneous interaction term, 

A a (x)A b (y)) = G(x-y)6 ab 6(x -yo) 
+ non- instantaneous, 

1 _o 1 



G(x-y)5° 



D[AJ • V D[A_l] • V 
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This expression describes an interaction generated by exchange of a Coulomb 

gluon dressed by the fluctuations of transverse gluon fields in the vacuum. 

By setting = we would return to the Laplace operator G = — 1/V 2 oc 
1/k 2 that corresponds, in the coordinate space, to the canonical Coulomb poten- 
tial 1/ |x — y|. For small vacuum fields, g s A±/X7 <C 1, expanding perturbatively 
the Coulomb propagator G in Eq. (jllj) to the second order in g s produces the 
one-loop anti-screening effect as stated in Eq. ()6bj) . 

If we take, however, gluon fields in the QCD vacuum as large as 

g s A ± /V ~ g s A ± -L~ 1 

(with L a spatial extent of the field), the perturbative expansion in <? S A_|_ of the 
denominators in Eq. (Jll[) is no longer justified. Under such circumstances a qual- 
itatively new phenomenon takes place namely, the Coulomb (ghost) propagator 
may become singular: 

(D[AJ-V)C = V 2 Q, + A ± ,VC ] = 0. (12) 

Appearance of a "zero-mode" solution Co in the external field is a sign of an 
infrared instability of the theory. 

An illuminating way to see the instability of the perturbative vacuum, and to 
re-derive the running of the coupling in QCD, is to look at the quantum correction 
to the vacuum energy Vq(H) = \H 2 in a constant chromo-magnetic field H. 

The one-loop corrected energy density reads 

Re V(H) = \H< + ig ,H)^ (h, Sf - i) * • V (H), (13) 

where // is the renormalization scale of the bare coupling, g s = g s ([i 2 ), and 
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b = lliV c /3 (gluodynamics). For relatively small fields, 



cr rj ( !6^ 2 

H < H ~ — p— exp 



5s(^ 2 ) V b 9 2 (fJ> 2 )J ' 
the potential Eq. 1)13(1 formally turns negative and develops a minimum — "true 
vacuum" (?) - corresponding to non-zero expectation value of the chromo- 
magnetic field. It was soon realized, however, that this new "vacuum" is unsta- 
ble fifty. 

In pure gluodynamics the vacuum correction can be computed by summing 
over the Landau levels of gluons moving in the external field: 



it p OO , 

6V{H) = 9s dp z J2Yl \l^sH{n + l-s z )+pl 

J n=0s z =±l 



(14) 



The real part of (|14[) yields l)13j). Thus, the asymptotic freedom in this approach 
can be seen as arising from the paramagnetic response of the vacuum. 

At the same time, we readily see that one specific state in Eq. (jl4jl namely, 
n = with gluon spin parallel to the direction of the field, s z = 1, gives an 
imaginary contribution: 

„ U f3sH 2tt2 2 

ImSV(H) = ^ J ^ dp z Im - g s H - iO = = -g- • V (H). (15) 

Non-zero imaginary part of the effective potential signals that the vacuum is 
unstable (decays with time) and does not correspond to the true ground state of 
the theory. Note also that n = corresponds to the Landau level of the largest 
transverse radius, i.e. to the infrared region, as had to be expected. 

Once again, we come to the conclusion that naive perturbative treatment of 
non-Abelian gauge fields is flawed due to problems in the infrared. 
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3.3 Gribov Copies 

In (124[) Gribov localized the problem. He showed that the three-dimensional 
transversality condition Eq. (JJJ) actually does not solve the problem of gauge 
fixing. 

Consider a group f2, and let oj{x) be a function with values in this group. The 
basic principle that defines the corresponding gauge theory is that the vector 
potentials A^(x) and 

AM(x) = u(x)A^x)uj~ 1 {x) + d^co(x)u;~ l (x) (16) 

describe physically identical fields. Therefore, to avoid multiple counting, one 
has to impose a gauge-fixing condition of the general form 

F(A^;x) = 0. (17) 

To do the job, Eq. (|17|) should have a unique solution for arbitrary A. This can 
be achieved for a variety of gauge fixings within perturbation theory, when the 
fields are weak. In general, however, solutions may appears copious since the 
gauge fixing condition Eq. (fTTj) constitutes a system of nonlinear equations for 
the function lo(x). 

Indeed, Gribov found that due to essential non-linearity of the gauge trans- 
formation Eq. (|16|). a "transverse" field potential satisfying Eq. Q may actually 
happen to be a pure gauge field which should not be separately counted as an 
additional physical degree of freedom. He explicitly constructed such "transverse 
gauge fields" for the SU(2) gauge group and showed that the uncertainty in gauge 
fixing arises when the effective magnitude of the field becomes large, 

A± ■ L ~ — , 

9s 
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or, in other words, when the effective interaction strength (QCD coupling) be- 
comes of the order of unity, that is, in the non-perturbative region. More precisely, 
he found that it happens exactly when the Faddeev-Popov operator acquires a 
zero mode solution Eq. Q12|) that is, as we discussed above, in the infrared region 
where the vacuum enhancement of the dressed Coulomb gluon propagator 
becomes catastrophically large. 

Thus, the "surface" (D[Aj_] • V)Co = in the functional Ax-space marks 
the border (the "Gribov horizon" ) beyond which the solution of the gauge-fixing 
equation Eq. (J7|) becomes copious. From this point of view, the fact that the 
Coulomb propagator develops singularity does not necessarily mean that the 
Faddeev-Popov ghost "rises from the dead" by pretending to propagate as a 
particle. It rather tells us that we have failed to formulate the quantum theory of 
interacting non-Abelian vector fields, to properly fix physical degrees of freedom. 

The existence of "Gribov copies" means that the standard Faddeev-Popov 
prescription for quantizing non-Abelian gauge theories is, strictly speaking, in- 
complete and has to be modified. 

A number of studies explored further the emergence of Gribov copies. It was 
found, both analytically (36) and numerically ()37[) that Gribov copies can exist 
even inside the Gribov horizon, and that a more narrow "fundamental modular 
region" had to be defined to avoid the problem. Nevertheless, recently it has been 
shown ( 38 ) that the copies within the Gribov horizon actually do not contribute 
to any expectation values, and thus for all practical purposes the original Gribov's 
recipe of constraining Faddeev-Popov determinant to positive values is correct. 
A promising attempt to implement the Gribov fundamental domain in the 5- 
dimensional formulation of gluodynamics can be found in (39). 
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3.4 Coulomb Confinement 

Gribov himself addressed the quest of possible modification of the QCD quanti- 
zation procedure in the second lecture "Quantization of non-Abelian gauge theo- 
ries" (|25j) . The paper (|4U|) under the same title based on the two Winter School 
lectures is now a universally accepted (though disturbing) truth and during 25 
odd years since its appearance in 1978 was being cited more than 660 times, with 
increasing frequency. 

To properly formulate non-Abelian field dynamics, Gribov suggested to limit 
the integration over the fields in the functional integral to the so-called funda- 
mental domain, where the Faddeev-Popov determinant is strictly positive (the 
region in the functional space of transverse fields before the first zero mode 
Eq. (|T2l emerges). 

Gribov produced qualitative arguments in favour of the characteristic modifica- 
tion of the gluon propagator, due to the new restriction imposed on the functional 
integral. Effective suppression of large gluon field results, semi-quantitatively, in 
an infrared singular polarization operator II oc k~ 2 , 

2 

D-\k) = k 2 + U(k 2 ) ~ k 2 + ^ . (18a) 

The gluon Green function coincides with the perturbative one at large momenta 
(small distances), D(k) oc k~ 2 but vanishes off at k = 0, instead of having a pole 
corresponding to massless gluons: 

D ^ K FT^' (18b) 

The new non-perturbative parameter a 2 in Eq. Q18|) has dimension (and the mean- 
ing) of the familiar vacuum gluon condensate, a 2 ~ (a s (F^) 2 ), that emerged in 
the context of QCD sum rules (41). 
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Literally speaking, the ansatz Eq. (|18bj) cannot be correct since such a Green 

function would violate causality. 6 In reality, the gluon (as well as the quark) 

propagator should have a more sophisticated analytic structure with singularities 

on unphysical sheets, which would correspond, in the standard field-theoretical 

language, to unstable particles. 

At the same time, the modification of the Coulomb (ghost) propagator due to 

restricting the functional integral to the fundamental domain $25$ resulted in the 

singular small momentum behaviour 

1 a 



G(k) <x 



N c gl k*' 



corresponding to a linear increase of the interaction energy at large distances 
R = |x — y| between colour charges, V(R) oc oR. 

The idea of confinement emerging from dressed Coulomb exchange was fur- 
ther explored by Zwanziger (|42§ who has recently shown (j43j) that for the static 
interaction potential V(R) the following inequality holds: 

V(R) < V Covl (R). (19) 

This means that if confinement exists in pure gluodynamics, it should arise 

already at the dressed one-gluon exchange level - - "no confinement without 

Coulomb confinement" (|43|). The inequality Eq. (|19|) appears inevitable: an 

inclusion of "quantum" (gluons, quarks) degrees of freedom can only screen, that 

is, suppress, the classical (Coulomb) interaction, as we have discussed in the 

beginning of this Section. 

Recent lattice studies of the correlator of timelike link variables in Coulomb 

gauge show that the Coulomb interaction energy of static quarks indeed 
6 It is unfortunate, therefore, that the form Eq. <18bt which Gribov suggested and discussed 

for illustrative purposes only, is often referred to in the literature as "the Gribov propagator" . 
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grows linearly at large quark separations. This simulation was designed to mea- 
sure the pure Coulomb energy, without admixture of additional "constituent" 
gluons, 

= g°(x) g a (y)* - (20) 

The corresponding string tension a turned out to be several times greater than 
the generally accepted value that originated from lattice measurements of the full 
static interaction energy. 

This fact is remarkable. It shows that the Coulomb exchange excites transverse 
vacuum gluons inclusion of which results in the energy gain. These gluons "shake" 
the initial string Eq. (|2U|) and soften it. 

In this picture the typical transverse size of the string grows logarithmically 
with quark separation (the so-called "roughening"). It is worthwhile to notice 
that this phenomenon is similar to the increase of the size of a hadron at high en- 
ergies as described by "Gribov diffusion" , multiplication of the number of gluons 
being the physical reason for both. 

3.5 Perturbative (Gluon) Confinement 

Motivated by Gribov ideas, J. Greensite and C. Thorn recently suggested and 
developed an interesting model that aims at reproducing the main features of the 
Coulomb confinement essentially perturbatively (45). 

They proposed the model in which gluons are arranged in chains that stretch 
between external charges. The Fock state wave function of such a chain looks 
like 

* chain [A} = <f^x.) A a ^( Xl ) A a ^( X2 ) ... A^+Hxn) q aN+1 (y)*o[A]. (21) 
The spatial distribution of the gluons is then determined by minimizing the total 
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energy. 

Qualitative explanation of the Coulomb confinement in terms of "the gluon 
chain model" of Greensite and Thorn runs as follows. In 't Hooft's large-iV c limit 
(N c —* oo, glN c fixed), the dominant diagrams are the planar ones in which each 
of the gluons interacts only with its nearest neighbors. (Note that the nearest- 
neighbor interaction can be purely Coulombic, i.e. perturbative.) This state 
can be approximately represented (modulo corrections) as a superposition 

of (N + 1) fundamental colour dipoles. As the distance between the quarks 
increases, so does the number of gluons (dipoles) N between them, so that the 
number of gluons per unit length turns out to stay finite, N/R = 1/1 =const, as 
Monte Carlo simulations showed. If -Egiuon is the kinetic plus nearest-neighbor 
interaction energy of each gluon, then the total energy of the system is 

E(R) = NE glnon = R = aR, (22) 

where the the role of the string tension is played by gluon energy per unit length, 
c = E g \ uon /t. 

It has been argued (j45j) that this model of confinement explains the Casimir 
scaling of the string tension observed on the lattice (the proportionality of the 
confining force between static sources in representation r of the gauge group to 
the quadratic Casimir operator C r of the representation) and yields the expected 
"center dependence" (the dependence of the string tension on the transformation 
properties of the center subgroup of the gauge group Q — "iV-ality"). 

Thus the linear growth of the interaction potential with distance is the con- 
sequence of the linear increase in the number of gluons excited by the static 
quark-antiquark interaction (|45[). It is important to stress that N linearly in- 
creasing with distance in the static problem, translates into the uniform rapidity 
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distribution of excited gluons in the relativistic situation when external sources 

fly away with the speed of light as, for example, in electron-positron annihilation 

e + e~ — > qq — > hadrons. 

Such a picture therefore directly relates the notion of a nearly constant confin- 
ing inter-quark force that is used to describe mass spectra of heavy quark bound 
states (|46j) . with well developed "string" phenomenology of multiparticle produc- 
tion in high energy interactions Q47|) . The latter successfully describes the gross 
features of both soft and hard phenomena. On the soft side we have "multiperiph- 
eral" hadroproduction in hadron-hadron reactions (Gribov plateau = Pomeron); 
on the hard side — the famous Feynman plateau, xD(x) ~const., characteris- 
tic for parton distributions that explain spacelike deep inelastic scattering (and 
particle content of timelike jets) at small Bjorken (Feynman) x. 

This important link is not a monopoly of the gluon chain model. Indeed, we 
seem to be retelling the good old Kogut-Susskind confinement story: gluon field 
lines forming a string, digging up quark-antiquark pairs from the vacuum, a'la 
Schwinger tunneling mechanism in a constant field, assembling colourless hadrons 
in the final state. There is a crucial difference, however. We were used to look at 
the Kogut-Susskind scenario as calling for essentially non-perturbative dynam- 
ics: strong fields, "superconductivity" vortex picture to explain one-dimensional 
nature of the "string", etc. Now we are lead to think that the same physical 
picture can be achieved by purely perturbative means, that is, by employing the 
language of the fundamental degrees of freedom of the theory — gluons — which 
gluons interact in a Coulomb-like (again, perturbative) fashion. 

't Hooft (48) has advanced the Greensite-Thorn arguments, and pointed out 
that the "gluon chains" cannot provide a complete set of states to confine the 
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sources. The problem becomes apparent if we turn to the real world with light 

quarks present. Indeed, by pulling apart two heavy quarks we expect to find 
ourselves holding two colourless D-mesons, in the end of the day. Obviously, 
this cannot be achieved without including light vacuum quarks in the game. 
(Formally speaking, the adjoint representation gluons cannot fully blanch colour 
charges in the fundamental representation — the external quarks.) Nevertheless, 
since gluon chain states appear to be energetically profitable, as recent lattice 
results have indicated (|45|). this model can serve as a good starting point, as a 
"first approximation" to the problem of quantitatively approaching the physics of 
hadrons. In particular, it may help to understand puzzling softness of the trans- 
formation of partons into hadrons. These phenomena that one observed studying 
energy and angular distributions of soft particles produced in hard interactions 
known as "local parton-hadron duality" (|49j) . for recent reviews see (|50U51|) . 

't Hooft's (unpublished) remark (|48|) bears an elegant strikingly simple title 
Perturbative Confinement. Looking back into historical perspective, we find this 
rather ironic, since this very endeavour — perturbative approach to the confine- 
ment problem — motivated the NPB referee to initially reject 7 the pioneering 
Gribov paper (|4U|) . 

Instability of the perturbative vacuum tells us that in order to have a chance 

to approach the true groung state of the theory by adiabatically switching on the 

interaction, one has to start from an excited state, in terms of non-interacting 

perturbative vacuum. Implementation of this idea is being developed by P. Hoyer 

and S. Peigne, who are trying to model non-trivial structure of the vacuum by 
7 The second reason being, that "the confinement problem had been already solved and isn't 

worth talking about." 
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explicitly adding condensate terms to perturbative quark and gluon propaga- 
tors (EH). 



3.6 OPE and Linear Potential 

Another way of approaching the problem of confinement is to start at short 
distances, use the operator product expansion (OPE) and to study the structure 
of power corrections. For the case of potential acting between massive quarks 
inside a colourless bound state, the expansion takes the form see for a 
recent survey, 

V(R) ~ \ 1 + 2_^a n ct s + c 3 - -if I, (23) 

where the second term in the brackets describes the higher order perturbative cor- 
rections, and the third one is the leading power correction. The latter is motivated 
by the OPE that predicts the leading non-perturbative power correction to the 
vacuum energy to scale in accordance with the dimension of the vacuum expec- 
tation value of the first gauge invariant vacuum operator, [(a s (F® u ) 2 }] = [-R -4 ]. 
Within perturbation theory, its origin can be traced to "infrared renormalons" , 
for a review see (55). 

What is most remarkable about the expression Eq. (|23|) in the context of our 
discussion is that it does not contain a term linear in R and thus contradicts the 
expected form 

V(R)~~ + *R, (24) 

borne out both by phenomenology of heavy quarkonia and by the lattice studies, 
as we have discussed above. We therefore come to the troubling conclusion that 
the OPE is inconsistent with the linear confining potential. 
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Yet another unexpected blow to the OPE ideology came from the direct state- 
of-the-art Monte Carlo simulation of the vacuum condensate in lattice gauge 
theory (|56[) which showed that the leading non-perturbative correction to the 
Wilson loop plaquette expectation value scales with lattice spacing as a 2 , instead 
of the OPE-blessed a 4 . 

Various options were suggested to avoid this contradiction: (see (|51)l and ref- 
erences therein). One can invoke new, genuinely non-perturbative, degrees of 
freedom as means for constructing missing dimension two vacuum condensates. 
The vocabulary of such approaches includes the notions of Dirac strings, (clus- 
tering) Monopoles, (percolating) Vortices and alike Q57[) . 

Alternatively, one may look at the QCD coupling itself; in other words, to 
search for 1/Q 2 corrections emerging from the dressed gluon propagator. 

3.7 Infrared-Finite QCD Coupling 

It is interesting that such terms naturally appear when one tries to implement 
the idea of analyticity in the Q 2 variable by assuming the existence of spectral 
representation for the running coupling (58,59). For example, a l/Q 2 -suppressed 
correction emerges as a result of simple removal of the Landau pole from the 
one- loop coupling (60): 



The ansatz Eq. 1)251) may turn out to be too simplistic. 8 Nevertheless, the idea of 
enforcing analyticity (causality) on a s turned out to be very efficient, see in partic- 
ular H61|) . The corresponding technology that improves perturbative expansions 

is known under the name of "Analytic Perturbation Theory" ()62|) . Impressive 
8 The sign of the 1/Q 2 term it produces does not match that of the lattice result 15411 . 




(25) 
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recent results of refurbishing of perturbative series can be found in (63). 

A QCD couling that stays finite all over the complex Q 2 -plane (as an analytic 
a s (Q 2 ) obviously does) allows one to peek into the strong interaction domain. 
One can aim at quantifying genuine confinement effects in QCD observables by 
linking power-suppressed non-perturbative contributions with momentum inte- 
grals (moments) of the coupling in the infrared region ( 59 ) . 

Markedly, the average value of the coupling that emerges from the study of the 
leading 1/Q power corrections to various jet shape observables in e + e~ annihila- 
tion and DIS turns out to be 

i r2 GeV 

ao = <a s (Q 2 )> = ^ F y o dQa s (Q 2 ) = 0.47 ±0.07. (26) 

(For a recent review see (|51|) and references therein.) Let us mention, in passing, 
that this number turns out to be suspiciously close to the corresponding integral 
over Eq. (|25|) . Even more interesting a "coincidence" is that the "measured" value 
Eq. (|26j) is comfortably above the so-called critical value of the strong coupling 
which is necessary, as we shall discuss in the next Section, to trigger the light 
quark confinement, according to Gribov, see Eq. (|28|) below. 

The quest for defining a s at large distances has a long and turbulent history, 
for reviews see (64,65,66). On one hand, it goes without saying that the "Landau 
pole" is an artifact. On the other hand, in QCD we don't have means for defining 
the true "physical coupling" unambiguously, in contrast with QED where a(0) is 
directly accessible via small angle scattering, Josephson effect and alike. On the 
phenomenological side, Mattingly and Stevenson (64) have assembled an impres- 
sive list of practical applications which consistently pointed at a s j-K = 0.2 — 0.3 as 
a reasonable magnitude of the "long-range" QCD interaction strength. The ap- 
plications they considered ranged from rather naive estimates of hadron-hadron 
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cross sections and form factors to the well elaborated Godfrey-Isgur relativized 

QCD quarkonium model that described quite successfully particle spectroscopy 
from pions all the way up to the T family. On the theory side, this quest as- 
cends to the notion of "effective charges" introduced by G. Grunberg back in 
1984 (67). In the context of the hunt for confinement effects, the concept of 
an infrared-finite coupling was suggested in (|68j) for the purpose of quantifying 
non-perturbative Aqqd/Mq terms in fragmentation functions of heavy quarks. 
Perturbative ideology behind this concept was laid down in (69). 

The model coupling Eq. (|^)) freezes at a constant value (a s (0) = An/bo). 
Though such a regime is often advocated in the literature, 9 we'd rather it vanished 
at the origin: any (even very weak) singularity at k 2 = of the dressed gluon 
propagator D(k) oc a s (k 2 )/k 2 would introduce unwanted long-range Van-der- 
Waals forces between hadrons. 

4 GRIBOV CONFINEMENT 

Returning to the task of constructing consistent QFT dynamics of non-Abelian 
gauge fields, we must conclude that, in spite of many attempts, the problem 
of Gribov copies ("Gribov horizon", "Gribov uncertainties") remains essentially 
open today. 

By the mid-80's Gribov decided, however, to change direction and not pursue 
pure gluodynamics. He did it not because of severe difficulties in describing the 
fundamental domain in the functional space: he always had his ways around tech- 
nical obstacles. Gribov convinced himself (though not yet the physics community 

at large) that the solution to the confinement problem lies not in the understand- 
9 for a recent emotional review see <!66t 
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ing of the interaction of "large gluon fields" but instead in the understanding 

of how the QCD dynamics can be arranged as to prevent the non-Abelian fields 

from growing real big. There was a deep reason for this turn, which he formulated 

in the following words: 

"I found I don't know how to bind massless bosons" 

(read: how to dynamically construct glueballs). 

As for fermions, there is a corresponding mechanism provided by the Fermi- 
Dirac statistics and the concept of the "Dirac sea". Spin-i particles, even massless 
which are difficult to localize, can be held together simply by the fact that, if 
pulled apart, they would correspond to the free-fermion states that are occupied 
as belonging to the Dirac sea. 

4.1 Perturbative Light Quark Confinement 

As a result of the search for a possible solution of the confinement puzzle Gribov 
formulated for himself the key ingredients of the problem and, correspondingly, 
the lines to approach it: 

• The question of interest is not of "a" confinement, but that of "the" con- 
finement in the real world, namely, in the world with two very light quarks 
- u and d - whose Compton wave lengths are much larger than the char- 
acteristic confinement scale (m q ~ 5 — 10 MeV <C 1 GeV). 

• No mechanism for binding massless bosons (gluons) seems to exist in Quan- 
tum Field Theory (QFT), while the Pauli exclusion principle may provide 
means for binding together massless fermions (light quarks). 

• The problem of ultraviolet regularization may be more than a technical trick 
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in a QFT with apparently infrared-unstable dynamics: the ultraviolet and 

infrared regimes of the theory may be closely linked. Example: the pion 

field as a Goldsone boson emerging due to spontaneous chiral symmetry 

breaking (short distances) and as a quark bound state (large distances). 

• The Feynman diagram technique has to be reconsidered in QCD if one 
goes beyond trivial perturbative correction effects. Feynman's famous ie 
prescription was designed for (and is applicable only to) the theories with 
stable perturbative vacua. To understand and describe a physical process in 
a confining theory, it is necessary to take into consideration the response of 
the vacuum, which leads to essential modifications of the quark and gluon 
Green functions. 10 

Existence of light quarks is crucial for the Gribov confinement scenario. It is 
clear without going into much mathematics that the presence of light quarks is 
sufficient for preventing the colour forces from growing real big: light quarks in 
the vacuum are eager to screen any separating colour charges and turn dragged 
apart heavy quarks into a pair of blanched Z?-mesons. 

The question becomes quantitative: how strong is strong? How much of a 
tension does one need to break the vacuum and organize such a screening? 

4.2 Supercritical Binding 

In a pure perturbative (non-interacting) picture, the empty fermion states have 

positive energies, while the negative- energy states are all filled. With account 

of interaction the situation may change, provided two positive- energy fermions 
10 The proper technology lies in a generalisation of the Keldysh diagram technique designed 

to describe kinetics out of equilibrium. 
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(quarks) were tempted to form a bound state with a negative total energy. In 
such a case, the true vacuum of the theory would contain positive kinetic energy 
quarks hidden inside the negative energy pairs, thus preventing positive-energy 
quarks from flying free. 

A similar physical phenomenon is known in QED under the name of super- 
critical binding in ultra- heavy nuclei. Dirac energy levels of an electron in an 
external static field created by the large point-like electric charge Z > 137 be- 
come complex. This means instability. Classically, the electron "falls onto the 
centre". Quantum-mechanically, it also "falls", but into the Dirac sea. 

In QFT the instability develops when the energy e of an empty atomic electron 
level drops, with increase of Z, below — m e c 2 . An e + e~ pair pops up from the vac- 
uum, with the vacuum electron occupying the level: the super-critically charged 
ion decays into an "atom" (the ion with the smaller positive charge, Z — 1) and 
a real positron 



Thus, the ion becomes unstable and gets rid of an excessive electric charge by 
emitting a positron ((7UJ). In the QCD context, the increase of the running quark- 
gluon coupling at large distances replaces the large Z of the QED problem. 

Gribov generalised the problem of supercritical binding in the field of an in- 
finitely heavy source to the case of two massless fermions interacting via Coulomb- 
like exchange. He found that in this case the supercritical phenomenon develops 
much earlier. Namely, a pair of light fermions interacting in a Coulomb-like 
manner develops supercritical behaviour if the coupling hits a definite critical 
value dZTJ) 



A z 



A Z -i + e 



+ 



for Z > Z crit . 




(27) 
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In QCD one has to account for the colour Casimir operator. Then the value of 
the coupling above which restructuring of the perturbative vacuum leads to chiral 
symmetry breaking and, likely, to confinement (see Q72|) and references therein), 
translates into 



This number, apart from being easy to memorize, has another important quality: 
it is numerically small. Gribov's ideas, being understood and pursued, offer an 
intriguing possibly to address all the diversity and complexity of the hadron world 
from within the field theory with a reasonably small effective interaction strength. 

4.3 Gribov Equation 

Gribov constructed the equation for the quark Green function as an approxima- 
tion to the general corresponding Schwinger-Dyson equation. This approxima- 
tion took into account the most singular (logarithmically enhanced) infrared and 
ultraviolet effects due to quark-gluon interactions and resulted in a non-linear 
differential equation which possesses a rich non-perturbative structure. 

An amazing simplicity of the Gribov construction makes one wonder, why such 
an equation had not been discovered some 20 years earlier when a lot of effort 
was applied in a search for non-perturbative phenomena of the superconductivity 
type in QED (Nambu-Jona-Lasinio; Baker-Johnson; Fomin-Miransky et al.). 

Take the first order self-energy diagram Ei(g): a fermion (quark/electron) with 
momentum q virtually decays into a quark (electron) with momentum q' and a 
massless vector boson (gluon/photon) with momentum k = q — q': 




N*-l 
2N C 



(28) 





k 2 + ie 



1 



(29) 
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with G and D the fermion and boson propagators, respectively. The correspond- 
ing Feynman integral diverges linearly at q' —* do. To kill the ultraviolet diver- 
gences (both linear and logarithmic), it suffices to differentiate it twice over the 
external momentum q. 

The first Gribov's observation was that 1/k 2 of the boson propagator happens 
to be the Green function of the four-dimensional Laplace operator, 

1 d 

d h A2 ; ■ = ~ 47r2 ^(g - g')» ^ = — , 

^{q-q') z + ie p dq^ 

where denotes the momentum differentiation. Therefore, the operation <9 2 
applied to Eq. (|29j) takes away the internal momentum integration and leads to 
an algebraic expression which is local in the momentum space, k = q — q' = 0, 



- for QED, 
* (30) 

C F — for QCD. 

This is the "Born approximation" . With account of higher order radiative correc- 
tions, the first thing that happens is that the bare fermion propagator Go dresses 
up, Go(q) — ► G(q), and so do the Born vertices 7 M — ► T^(q,q,0). The second 
crucial observation was that the exact vertex function T^q^q — k,k) describing 
the emission of a zero momentum vector boson, fc„ = 0, is not an independent 
entity but is related with the fermion propagator by the Ward identity, 

r M (g,g,0) = -d II G-\q). (31) 

This statement is literally true in Abelian theory (QED), and, after some reflec- 
tion, can be made true in the non-Abelian case (QCD) as well. 

Thus, we arrived to the Gribov equation for the quark Green function H71II72|) 

dlG-\q)=gd li G-\q)G(q)d ll G-\q) + ... , (32) 

where . . . stand for less singular O (g 2 ) integral terms. 
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In principle, the whole PT series expansion may be constructed for the right 

hand side in terms of exact Green functions (and their momentum derivatives). 

In particular, with account of the first subleading terms Eq. l|3~2|) becomes an 

integro-differential equation and its r.h.s. may be represented in the following 

graphic form: 



dlG-\q) 




+ 



with black dots standing for the momentum gradient of the inverse propagator - 
the exact zero- momentum boson emission vertex ( 69 ) . Yet another set of higher 
order corrections makes the coupling run, g — > g(q 2 )- In the \q 2 \ — > oo limit the 
QCD coupling vanishes due to the asymptotic freedom, and Eq. (|32j) turns into 
the free equation, d^G^ 1 = 0, whose general solution has the form 



G-\q) = Z Q l 



4 A 



(m - q) + — + — 



(q = 1^) • (33) 



This general perturbative solution has two new arbitrary parameters v\ and V2 in 
addition to the familiar two (bare mass mo and the wave function renormaliza- 
tion constant Zq), since the master equation is now the second order differential 
equation, unlike in the standard renormalization group (RG) approach. 

The new terms are singular at q 2 — > 0. Therefore in QED, for example, we 
simply drop them as unwanted, thus returning to the RG structure. Such a 
prescription, however, exploits the knowledge that nothing dramatic happens 
in the infrared domain, so that the real electron in the physical spectrum of the 
theory, whose propagation we seek to describe, is inherently that very same object 
that we put into the Lagrangian as a fundamental bare field. 

Not so clear in an infrared unstable theory (QCD). Here we better keep all four 
terms in Eq. (|3l?|). wait and see. 
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At large virtualities, \q 2 \ S> m?, the two additional terms can be looked upon 

as power suppressed corrections that emerge due to non-trivial structure of the 

QCD vacuum. The corresponding dimensional parameters can be directly linked 

with the non-perturbative vacuum condensates, 

which constituted the core of the "ITEP sum rules" famous for successful phe- 
nomenology of low- lying hadron resonances, see 1)73(1 . 

Moreover, in the finite momentum region, \q 2 \ ~ m 2 , where all four terms have 
to be treated on the same footing, Gribov found bifurcation - a non-perturbative 
solution - emerging in Eq. ()32j) if the coupling in the infrared region exceeded 
the critical value Eq. (|28|) . The new phase corresponds to spontaneously broken 
chiral symmetry. This means that given a supercritical coupling in the infrared, 
the quark Green function may possess a non-trivial mass operator even in the 
chiral limit of vanishingly small bare (ultraviolet) quark mass tuq — > (|711I72|) . 

4.4 Critical Coupling and Chiral Symmetry Breaking 

In this section we suggest a simple algebraic exercise that should help a reader 
to see the origin of the critical coupling and to appreciate the Gribov equation 
as a tool for grasping dynamical breaking of the chiral symmetry. 
Let us introduce 

A u (q) = d u G-\q)-G{q). (34) 

Observing that 

dlG' 1 - d v {dyG' 1 ■ G) = -d v G- x ■ d v G = (d^G' 1 ■ G) (d v G- 1 ■ G) , 
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the master equation Eq. (|32j) can be cast in the following compact form, 



d v A v + {l-g) A V A U = 0. 



(35) 



Using the standard representation for the fermion propagator in terms of the 
wave function renormalization factor Z and the running mass m, 

G-\q) = Z-\0 [m(0 -q], £ = lng^ In v^, 



and introducing anomalous dimensions 



Z 



m 
m 



f = 



(36) 



it is straightforward to derive an explicit expression for the "vector field" A u : 



A, 



m 



Here 



o-y^ + (1 - T) n v (r m - T)hn u 



1 



1 + hm/q 
1 — (m/q) 2 



(37) 



n u = — , n =1; a V{i = -icivly. ~ lulu), cr^a^ = -3g up . 
q z 

In the region of relatively large momenta, q 3> m, we have 

TTl f TYl 2 \ 

q ■ A v = [ o^n* + (1 - T)n v ] + — ( 7 „ - T m hn v ) + O 
Applying this approximation to Eq. (|35j) gives 



q 2 d v A v ~ 2(i-r)-r 



Em ~\~ E m + 2 



/n 



• n; 



(l-E) 2 -3 + 2(i-r)(i-r r 



m 



• n. 



(38a) 
(38b) 



Equating the scalar and h terms in Eq. (|35|) produces then a system of coupled 
differential equations for anomalous dimensions: 



f = (l- 9 )(l-E) 2 + 2(l-E)-3(l-£), 



(39) 



(E 2 „ + 2)+2(l- 5 )(l-E)(l-r m ). 



(40) 
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It is its stable point T = T m = which determines the ultraviolet anomalous 

dimensions T* and T^. The first equation T = is self-contained and results in 

(i-,)(i-r*) = V3(i- 5 ) 2 + i-i, r* = 2 " VW-jp ±z9 t (41) 

where the sign of the square root was chosen such as to select among the two 
fixed points in Eq. (|39|) the stable one. 11 The anomalous dimension of the mass 
in Eq. (|4U|) is driven by that of the wave function, since from T m = we have 

r^ ± = -(i - g )(i - r) ± V[(i-s)(i-r) + i] 2 -3. (42) 

Substituting T from Eq. (|41|). for the stable point we finally obtain 

= = 1 _V3(l- 5 ) 2 + l + \/3(l-5) 2 -2. (43) 

Expressions Eqs. H41|) - (|43|) ar e non-perturbative since they have emerged from 
non-linear Eq. (|32j) . The first term of the series expansion in g, 

T*(g) = l -g + (g 2 ) , T* m (g) =-\g + (g 2 ) , 

coincides with the known one-loop perturbative expression for the anomalous 
dimension of the (Feynman gauge) fermion wave function and of the (gauge in- 
variant) mass operator, respectively. 

Now we can readily see how the critical coupling emerges that have been an- 
nounced above in Eq. (|27jl. Indeed, when g reaches the value g cr it = 1 — y^§, the 
running mass becomes complex due to the term y3(l — g) 2 — 2 in Eq. (|43[). 

This signals instability. In fact, at this point the two anomalous dimensions of 

Eq. (|4*2*)) become equally important, ReT^ n _ ! _(g = g cr it) = -(V3- !)• The mass 

operator m(£) remains real but ceases to be monotonic and starts to oscillate with 
n The second solution corresponds to F*(g = 0) —4 and describes renormalization of the 

power suppressed Z~V|/g 4 term of the free fermion propagator Eq. I)33|l . 
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£ = In q. It is this oscillatory behaviour that allows one to construct a symmetry 

breaking solution corresponding to mo = 0. Its mass operator is regular at q = 
and decays fast in the ultraviolet, m(£) oc exp(— 2£) oc l/q 2 , where the chiral 
symmetry gets dynamically restored. 

For explicit construction of this chiral symmetry breaking solution, its proper- 
ties and physics of accompanying Goldstone states see (J72J). 

4.5 Quarks, Pions and Confinement 

As far as confinement is concerned, the approximation Eq. (|32|) turned out to be 
insufficient. A numerical study of the Gribov equation carried out by C. Ewerz 
showed (|74[) that the corresponding quark Green function does not possess an 
analytic structure that would correspond to a confined object. 

Given the dynamical chiral symmetry breaking, however, the Goldstone phe- 
nomenon takes place bringing pions to life. In his last paper (75 ) Gribov argued 
that the effects that Goldstone pions induce, in turn, on the propagation of quarks 
is likely to lead to confinement of light quarks and, as a result, to confinement of 
any colour states. 

The approximate equation for the Green function of a massless quark, which 
accommodates a feed-back from Goldstone pions reads (|75|) 

dlG~\q) = g(q)d IM G~ 1 (q)G(q)d^ 1 (q) 

-leJrp {^G-Hq)} G(q) {i^G-\q)}. (44) 

It is important to notice that since pions have emerged dynamically in the the- 
ory, their coupling to quarks is not arbitrary but is tightly linked with the quark 
propagator itself (search for an anti-commutator of 75 with G~ l in Eq. (144(1 ). 
Moreover, the pion-axial current transition constant f v is not arbitrary either, 
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but has to satisfy a definite relation which, once again, is driven by the behaviour 
of the exact quark Green function: 



The second of the two papers (|72U75|) concluding Gribov's study of the light- 
quark supercritical confinement theory remained unfinished. It ends abruptly in 
the middle of the discussion of the most intriguing question, namely, what is the 
meaning, and practical realization, of unitarity in a confining theory. 

The modified Gribov equation Eq. Q44|) still awaits a detailed study aiming at 
the analytic structure of its solutions. 

4.6 Gluon Sector 

Another important open problem is to construct and to analyse an equation for 
the gluon similar to that for the quark Green function, from which a consis- 
tent picture of the coupling g{q) rising above the critical value in the infrared 
momentum region should emerge. 

The difficulty with the gluon sector of the theory lies in separating the running 
coupling effects from an unphysical gauge dependent phase that are both present 
in the gluon Green function. To analyse renormalization of the gluon Gribov has 
used in (|72[) the Duffin-Kemmer formalism (also known as "linear formalism") 
which treats the gluon potential A® and the field strength as independent 
variables. Renormalization properties of the gauge invariant correlator {FF) 
gives then a direct access to the running coupling. The Duffin-Kemmer technique 
being plagued with artificial divergencies, this attempt did not result however in 
an equation for ct s in a closed form. 



J 7T 




(45) 
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The solution may lie in using the "background gauge" or in exploiting the uni- 

tarity motivated "pinch" technique developed by J. Papavassiliou, N.J. Watson 

and others (see (|75)l and references therein), as in both approaches the QCD 

coupling is directly related with the renormalization of the gluon Green function. 

5 QED AT SHORT DISTANCES 

At the same time, to construct a smart equation for the running coupling poses 
no problem in an Abelian gauge theory where the boson propagator is gauge 
invariant. 

Gribov has carried out this programme for QED. Here the coupling becomes 
supercritical at extremely short distances. Gribov's analysis aimed at resolving 
the long-standing problem of the "Landau pole" in the running QED coupling. 
The "Landau ghost" in the photon propagator at academically large momenta 
seems to be a formal problem so long as QED is actually a part of a broader field 
theoretical scheme. Nevertheless, its resolution is extremely important since it 
demonstrates what type of new, non-perturbative, phenomena in QFT one might 
expect when the strength of the coupling becomes large (dynamical symmetry 
breaking and appearance of Goldstone states, condensation, confinement, etc.). 

Gribov did not finish the paper entitled "Quantum Electrodynamics at Short 

Distances". Handwritten notes he left behind, deciphered and translated into 

English by J. Nyiri, appeared in 1)77(1 . It is important to mention that no attempt 

has been made by the editorial team to "debug" these notes. 12 

12 This was done on purpose. By correcting omnipresent mistakes/misprints in obvious places 

the editors could have misled a potential reader into putting too much trust into not-so-obvious 

derivations and conclusions. As the matter stands now, it should be not read but rather worked 

through with a pencil and an A4-pad at hand. We strongly encourage to do so everyone willing 
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Here we will reproduce essential steps of the part of Gribov's analysis that 

directly addressed the "Moscow zero" problem. By so doing we hope to be able 
to ignite enough curiosity so that our readers will be tempted to look them- 
selves into more delicate issues (such as appearance of light super-bound states, 
condensation, the Goldstone phenomenon) that were also considered in l|77|). 

5.1 Equation for Photon Polarization Operator 

Electron loop Feynman diagrams for the polarization operator H^^k) for a pho- 
ton with momentum k diverges quadratically. Therefore in order to obtain a 
convergent integral (and thus a finite answer) we need to differentiate it thrice 
over k. At the same time, conservation of current dictates 

n^(fe) = {k^K - g^k 2 ) ■ n(0, £ = ink, (46) 

where the factor II diverges logarithmically in the ultraviolet; in the Born ap- 
proximation, 

Ho = e\ ■ const • In ■ 



Auv 



k 

Let us apply successively two differential momentum operations to Eq. (|46jl : 



a 8 « 



=0„ 3fc„n(0 =3(4n+n) (47) 



(recall that the dot marks ^-derivative). To get rid of the remaining logarithmic 
divergence it suffices to differentiate Eq. (|4*7|l over £. This is equivalent to applying 
the operator k a d a which produces 



ll + I n = %r*-'V u V(*->- 



to play a i-of-a-Rosetta-Stone quest. 
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The l.h.s. can be expressed directly through the running QED coupling. Indeed, 

observing that, by definition, 

a(0 e 2 (e) e 2 1 



9(0 



vr " 4vr 2 4vr 2 1 + 11(f) ' 



we may write 




d 1 d 2 \ 1 (fc^ CT ) 2vr 2 

^ + 4^J- = — --^-^n^ib). (48) 

As for the r.h.s. of Eq. (|48|) . we treat IL^ as a sum of Feynman diagrams and 
apply the first two differentiations diagrammatically. This results in 

V . - , . 

= ^ - k -7 \ ;') (49) 

e .. - 



where crosses mark momentum derivative of the bare electron propagator, 

d„G (q + k) = <9 M ) ^ = Go^Go; x = -c^G 1 = 7m . (50) 

mo - (q + fe) 

(Integral over the loop momentum q diverges only logarithmically which justifies 
integration by parts used to derive Eq. (|49|).) After differentiation, the overlap- 
ping divergencies disappear that plagued the original Schwinger-Dyson equation 
for the polarization operator. Therefore, higher order corrections combine to 
renormalize four fermion propagators, Go(q) — > G(q), both external vertices, 
7 M — > T(q, q + k,k), as well as the crosses: 

x=^ = -^-G^(q) -> ~G-\q)=T^q,q,0) = ?. (51) 
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We thus arrive at the renormalized graph 



1 k 

— d^d u U^(k) = -- 

e o 
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+ A. 



(52) 



For the correction term A = O (g) a perturbative series expansion can be con- 
structed in terms of two-particle irreducible multi-loop graphs built of exact prop- 
agators and interaction vertices, 

q 

k 





+ -4 / V- +0{g 2 ). 



A 



(53) 

Feynman integral corresponding to the leading contribution in Eq. (|52[) reads 

where we have combined the adjacent dot-operator * = d fJ ,G^ 1 (q) and the Green 
function G(q) into of Eq. (J3U). 

We are interested solely in the logarithmically divergent contribution to Eq. (|54[) 
since finite pieces, after applying the last differentiation (k^d^) in Eq. (|48|). will 
produce but negligible corrections that are power suppressed for large external 
photon momenta, k S> m. Bearing this in mind, in the relevant integration region 
k <C q <C Auv we can omit k in the loop propagators, G(q + k) — G(q), and in 
the exact photon vertices, T^(q,q+k,k) ~ F^(q,q,0) = —d^G^ 1 (q), 
In the large momentum region we invoke Eq. (|38b|) to obtain 

4 



Tr 



AJq)AJq) A u {q)A u {q) 



([i-r] 2 -3)' 
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Transforming the integration phase space, 



(2ir)Hq 4 8vr 2 
we derive 



-^d (l d v n. lw {k) ~ / d£ 

°0 Jink 



4 

Substituting into Eq. (|48|) and taking the ^-derivative we finally arrive at 



)~ g * -^([l-r] 2 -3) 2 . (55) 



d 1 d 2 \ 1 I 



This second order differential equation generalises the standard RG equation, 

with the /^-function expressed via anomalous dimension of the (Feynman gauge) 
electron wave function, T = T*(g). 

5.2 Small and Large Coupling Regimes 

In the small coupling regime we have r = O (g) <C 1 and Eq. ((531) gives 

-i -i 2 

which is the standard one-loop expression for the running QED coupling that 
formally develops "Landau pole" in the deep ultraviolet, k ~ m e exp(4^ • 137). 
However, for g = O (1) there is no reason to neglect T. The latter is given by 

Y{g) = l 3(1 ~ 9) . (56) 

W l + Vl + 3(l- 5 ) 2 

From this representation, which is equivalent to the original Eq. (|41|) . we observe 
that r is monotonically increasing with g, crosses unity at g = 1 and tends to 

T(g) 9 =°° (1 + V3) - g-' + O (<T 2 ) . 
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We observe that in the large coupling limit the r.h.s. of Eq. (|55|) becomes small, 

(A. I i (2y / 3) 2 _2 

resulting in logarithmic increase of the coupling at large photon virtualities: 

g (£) ~ 2-£ = ln|/c 2 | » 1. (57) 

Landau singularity has disappeared. Gribov explained this phenomenon as a 
compensation between the contributions to vacuum polarization from magnetic 
moment of electron [(cr^ u ) 2 = —3] by that of its charge [(1 — T) 2 ] in the non- 
perturbative regime of large anomalous dimensions. 

Numerically large coupling j> 1 does not mean that the interaction becomes 
really strong. If this were the case, the approximation based on neglecting higher 
order terms A in Eq. Q52|) that was adopted in Eq. (|55[) would have been under- 
mined. Gribov argued that since adding a photon is accompanied by additional 
factors in each vertex, the relative magnitude of radiative corrections, in spite 
of g 3> 1, may turn out to be finite, g ■ (A^) 2 = O (1). 

In fact, there is no need to analyse the structure of A for g 3> 1, which regime is 
actually of little interest. Indeed, as have we discussed above, already at g = O (1) 
a new phenomenon occurs, which is supercritical binding of fermion pairs. It 
leads to appearance of scalar /pseudoscalar "mesons" and changes drastically the 
behaviour of the Abelian theory in the ultraviolet. 

5.3 Composite Higgs as a Super-critically Bound it Pair 

As an interesting byproduct of his supercritical quark confinement in the infrared, 
Gribov considered in 1994 an application of the supercritical picture to the weak 
sector of the Standard Model in the ultraviolet. 
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It is widely believed that the couplings of the basic interactions U (1)y> SU (2) wea k 

and 5'C/(3) s trong merge at some Grand Unification scale, above which the under- 
lying dynamics becomes essentially different. However, if there is a "Super-Great 
Desert" instead, the Abelian hyper-charge interaction constant gy keeps growing 
with scale leading to spontaneous symmetry breaking, non-zero vacuum expec- 
tation value V, generation of W/Z masses and appearance of Higgs boson as a 
super-bound state of top quarks. All these phenomena being inter-connected, this 
scenario allows one to relate m# and mt with the symmetry breaking parameter 
V ~ 246 GeV, the Weinberg angle and the values of the couplings a s and a e .m. 
at the top quark mass scale (J75J). 

The mass of such composite Higgs particle (whose binding is similar to that 
of a deuteron in the zero-radius limit of nuclear forces) would have been ran = 
(170±7) GeV, given the experimental top quark mass mt ~ 175 GeV (and varying 
a s between 2m t and \m t ). 

6 CONCLUSIONS 

Hadron phenomenology has accumulated a very impressive dossier of puzzles 
and hints, ranging from unexplained regularities in hadron spectroscopy to soft 
"forceless" hadroproduction in hard processes. 

The reason why one keeps talking, 30 years later, about puzzles & hints, about 
constructing QCD rather than routinely applying it, lies in the conceptually new 
problem one faces when dealing with a non- Abelian theory with unbroken sym- 
metry. We have to learn to master a Quantum Field Theory whose dynamics is 
intrinsically unstable in the infrared domain so that the objects belonging to the 
physical spectrum of the theory (supposedly, colourless hadrons, in the QCD con- 
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text) have no direct one-to-one correspondence with the fundamental fields the 

microscopic Lagrangian of the theory is made of (coloured quarks and gluons) . 

In these circumstances we don't even know how to formulate at the level of the 

microscopic fields the fundamental properties of the theory, such as conservation 

of probability (unitarity) and analyticity (causality). Indeed, 

• What does Unitarity imply for confined objects? 

• How does Causality restrict quark and gluon Green functions and their 
interaction amplitudes? 

• What is the Mass of an INFO - [well] Identified [but] Non-Flying Object? 13 

Understanding the confinement of colour remains an open issue. The very 
problem can be formulated in various terms, ranging from a 10 6 $ worth a rigorous 
mathematical proof of the existence of a mass gap in pure gluodynamics all 
the way to developing down-to-earth but theory motivated practical recipes for 
cooking the spectrum of hadrons and predicting their production properties. 

Among many a contribution to the problematics of non-Abelian gauge theories 
- proper quantization of Yang-Mills fields, the origin of asymptotic freedom, the 
nature of instantons and physics of quantum anomalies, etc., - the key ingredi- 
ent of the Gribov conception of Quantum Chromodynamics was setting up the 
problem of the confinement of colour as that of light quarks. 

Gribov works on gauge theories and, in particular, all his papers, talks and 

lectures devoted to anomalies and the QCD confinement (including the lectures 
13 The issue of quark masses is particularly damaging since a mismatch between quark and 

hadron thresholds (whatever the former means) significantly affects predicting the yield of heavy- 
flavoured hadrons in hadron collisions, markedly at LHC. 
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that were translated into English for the first time) were collected and recently 

published 14 in Gauge Theories and Quark Confinement Q79j) . 

From these papers, an interested reader will be able to follow the derivation 

of the Gribov equation and to study the properties of its (perturbative and non- 

perturbative) solutions, as well as to formulate and pursue the open problems 

awaiting analysis and resolution. Pedagogical lectures he gave in 1992 in Or- 

say (80) will give an opportunity to grasp the physical picture of the supercritical 

binding which includes an anti-intuitive notion of an "inversely populated" Dirac 

sea and to think about phenomenological aspects of the light quark confinement 

scenario. 

Needless to say, the Gribov approach and recent developments reviewed here 
did not yet provide a consistent recipe for dealing with QCD at large distances. 
Nevertheless, an impressive progress has been made, and Gribov's ideas continue 
to inspire many of those who have contributed to it. 

Literature Cited 

1. Gribov VN. Nucl. Phys. 22:249 (1961) 

2. Mandelstam S. Phys. Rev. 112:1344 (1958) 

3. Regge T. Nuovo Cimento 14:951 (1959); ibid., 18:947 (1960) 

4. Frautschi SC, Gell-Mann M, Zachariasen F. Phys. Rev. 126:2204 (1962) 

5. Gribov VN. Zh.E.T.F. 41:667 (1961) [Sou. Phys. JETP 14:478 (1962)] 

6. Froissart M. Phys. Rev. 123:1053 (1961) 

7. Gribov VN, Pomeranchuk IY. Phys. Rev. Lett. 8:343 (1962); Zh.E.T.F. 

42:1141 (1962) [Sou. Phys. JETP 15:788L (1962)] 
14 The book can be ordered at http://www.prospero.hu/gribov.html 



Gribov Conception of Quantum Chromodynamics 55 

8. Chew GF, Frautschi SC. Phys. Rev. Lett. 7:394 (1961) 

9. Gribov VN. Zh.E.T.F. 41:1962 (1961) [Sov. Phys. JETP 14:1395 (1962)] 

10. Gribov VN. Zh.E.T.F. 42:1260 (1962) [Sov. Phys. JETP 15:873 (1962)] 

11. Mandelstam S. Nuovo Cimento 30:1148 (1963) 

12. Gribov VN. Space-time description of the hadron interaction at high energies. 
In Proc. 1 st ITEP School, Vol. 1 Elementary particles, p. 65. Moscow: ITEP 
(1973). Published in English in (J73J), pp. 3-27, |hep-ph/0006158| 

13. Feynman R. What neutrino can tell about partons. In Conference Proc. 
Neutrino-72, Hungary, 1972, Vol. II, p. 75 

14. Gribov VN. Zh.E.T.F. 57:1306 (1969) [Sov. Phys. JETP 30:709 (1970)] 

15. Ioffe BL. Phys. Lett. B30:123 (1969) 

16. Abramovsky VA, Gribov VN, Kancheli OV. Yad. Fiz. 18:595 (1973) [Sov. J. 
Nucl. Phys. 18:308 (1974)] 

17. Gribov VN. Zh.E.T.F. 56:892 (1969) [Sov. Phys. JETP 29:483 (1969)] 

18. Gribov VN, Migdal AA. Zh.E.T.F. 55:4 (1968); Zh.E.T.F. 55:1498 (1968) 
[Sov. Phys. JETP 28:784 (1969)]. In Regge theory of low-px hadronic inter- 
actions, ed. Caneschi L, p. 90. North-Holland (1989), 410 pp. 

19. Gribov VN. Yad. Fiz. 17:603 (1973) [Sov. J. Nucl. Phys. 17:313 (1973)] 

20. Gribov VN. Nucl. Phys. B 106:189 (1976) 

21. Fadin VS, Kuraev EA, Lipatov LN. Phys. Lett. B60:50 (1975); Zh.E.T.F. 
72:377 (1977) [Sov. Phys. JETP 45:199 (1977)] 

22. Gribov VN. Kharkov Physico- Technical Institute preprint (1969, unpub- 
lished). Published in English: Theory of complex angular momenta. Gribov 
lectures on theoretical physics. Cambridge University Press (2002) 

23. Gribov VN, Lipatov LN, Frolov GV. Yad. Fiz. 12:994 (1970) [Sov. J. Nucl. 



56 Dokshitzer Sz Kharzeev 

Phys. 12:543 (1971)] 

24. Gribov VN. Instability of Non-Abelian Gauge Theories and Impossibility of 
Choice of Coulomb Gauge. In Proc. XII Winter School, Vol. 1 Elementary 
Particle Physics, p. 147. Leningrad: LNPI (1977). Published in English in 
(|73)) . pp. 258-270 

25. Gribov VN. Quantization of Non-Abelian Gauge Theories. In Proc. XII Win- 
ter School, Vol. 2 High Energy Physics, p. 64. Leningrad: LNPI (1977). Pub- 
lished in English in J7HJ), pp. 233-255 

26. Gross DJ, Wilczek F. Phys. Rev. Lett. 30:1343 (1973); Politzer HD. Phys. 
Rev. Lett. 30:1346 (1973) 

27. Khriplovich IB. Yad. Fiz. 10:409 (1969) [Sov. J. Nucl. Phys. 10:235 (1970)] 

28. Landau LD, Abrikosov A, Khalatnikov I. Nuovo Cimento Suppl. 3:80 (1956) 
In Collected Papers of L.D. Landau, ed. D ter Haar, p. 701. Pergamon Press 
(1965), 836 pp. 

29. Pomeranchuk IY. Akad. Nauk USSR 103:1005, 104:51, 105:461 (1955); 
Abrikosov AA, et al. Phys. Rev. 111:321 (1958) In Collected Papers of L.D. 
Landau, ibid., p. 761 

30. Kallen G, Pauli W. Math.-fys. Medd. 30: No. 7 (1955) 

31. Dyson FJ. Sci. Amer. 199:74 (1958) 

32. Landau LD. In Theoretical Physics in the Twentieth Century, a Memorial 
Volume to Wolfgang Pauli, ed. Fierz M, Weisskopf VF, p. 245. New York: 
Interscience (1960), 328 pp. In Collected Papers of L.D. Landau, ibid., p. 800 

33. Christ NH, Lee TD. Phys. Rev. D 22:939 (1980) 

34. Batalin IA, Matinyan SG, Savvidy GK. Sov. I. Nucl. Phys. 26:214 (1977); 
Matinyan SG, Savvidy GK. Nucl. Phys. B 134:539 (1978); Savvidy GK. Phys. 



Gribov Conception of Quantum Chromodynamics 57 
Lett. B 71:133 (1977) 

35. Nielsen NK, Olesen P. Nucl. Phys. B 144:376 (1978) 

36. Semenov-Tyan-Shanskii M, Franke V. Zapiski Nauch. Sem. Leningradskogo 
Otdeleniya Matematicheskogo Instituta im. V.A. Steklova AN SSSR 120:159 
(1982). English translation: Plenum Press (1986) 

37. Mandula JE, Ogilvie MC. Phys. Rev. D 41:2586 (1990); Marinari E, Parinello 
C, Ricci R. Nucl. Phys. B 362:487 (1991) 

38. Zwanziger D. Phys. Rev. D 69:016002 (2004) 

39. Baulieu L, Zwanziger D. Journal of High Energy Physics 0108:015 (2001), 
hep-th/0107074 



40. Gribov VN. Nucl. Phys. B 139:1 (1978) 

41. Shifman MA, Vainshtein AI, Zakharov VI. Nucl. Phys. B 147:385 (1979); for 
a recent review, see Ioffe BL. Phys. Atom. Nucl. 66:30 (2002) 

42. Zwanziger D. Nucl. Phys. B 518:237 (1998) 

43. Zwanziger D. Phys. Rev. Lett. 90:237 (2003) 

44. Cucchieri A, Zwanziger D. Nucl. Phys. Proc. Suppl. 119:727 (2003); Greensite 
J, Olejnik S. Phys. Rev. D 67:094503 (2003) 

45. Greensite J, Thorn CB. Journal of High Energy Physics 0202:014 (2002), and 
references therein. 

46. Bali GS. QCD Potentiology. In Quark confinement and the hadron spectrum, 
p.209. Wien (2000), |hep^ph/0010032| 

47. Andersson B, Gustafson G, Peterson C. Phys. Scripta 19:184 (1979); for a 
detailed review see Andersson B. The Lund Model. Cambridge Monogr. Part. 
Phys. Nucl. Phys. Cosmol. 7:1-471 (1997) 

48. 't Hooft G. Perturbative Confinement, hep-th/0207179 (2002) 



58 Dokshitzer <fe Kharzeev 

49. Dokshitzer YL, Troyan SI. Asymptotic Freedom and Local Parton-Hadron 

Duality in QCD Jet Physics. IN Proc. XIX Winter School, Vol. 1, p. 144 
(in Russian). Leningrad: LNPI (1984); Azimov YI, et al.Zeit. Phys. C 27:65 
(1985); Dokshitzer YL, Khoze VA, Mueller AH, Troyan SI. Basics of Per- 
turbative QCD , ed. J. Tran Thanh Van. Gif-sur-Yvette: Editions Frontieres 
(1991) 

50. Khoze VA, Ochs W. Int. J. Mod. Phys. A 12:2949 (1997) 

51. Dokshitzer YL. QCD Phenomenology. In Proc. European School of High- 
Energy Physics, Pylos, Greece, 25 Aug - 7 Sep 2002, hep-ph/ 0306287 

52. Hoyer P. |hep-ph/0203236| (2002); Hoyer P, Peigne S. |hep-ph/0304010| (2003) 

53. Balitsky II. 1983. Nucl. Phys. B 254:166; Aglietti U, Ligeti Z. Phys. Lett. 
B364:75 (1995); Akhoury R, Zakharov VI. Phys. Lett. B438:165 (1998) 

54. Zakharov VI. Nonperturbative Match of Ultraviolet Renormalon. 
hep-ph/0309178 (2003) 

55. Beneke M. Phys. Rep. 317:1 (1999) 

56. Burgio G, Di Renzo F, Marchesini G, Onofri E. Phys. Lett. B422:219 (1998), 
hep-ph/9706209 



57. Zakharov VI. Fine Tuning in Lattice SU(2) Gluodynamics versus Continuum 
Theory Constraints. In Int. Conf.: I.Ya. Pomeranchuk and Physics at the 
Turn of Centuries, Moscow, Russia, Jan 24-28, 2003, hep-ph/ 0306262 

58. Ball P, Beneke M, Braun VM. Nucl. Phys. B 452:563 (1995) 

59. Dokshitzer YL, Marchesini G, Webber BR. Nucl. Phys. B 469:93 (1996) 

60. Shirkov DV, Solovtsov IL. Phys. Rev. Lett. 79:1209 (1997); Phys. Part. Nucl. 
32SL48 (2001); Grunberg G. Journal of High Energy Physics 0108:019 (2001) 

61. Geshkenbein BV, Ioffe BL. JETP Lett. 70:161 (1999), |hep-ph/9906406l 



Gribov Conception of Quantum Chromodynamics 59 
Geshkenbein BV, Ioffe BL, Zyablyuk KN. Phys. Rev. D 64:093009 (2001), 

|hep-ph/0104048| 

62. Shirkov DV. Eur. Phys. J. C 22:331 (2001), hep -ph/0107282 



63. Howe DM, Maxwell CJ. All-Orders Infrared Freezing of Observables in Per- 
turbative QCD. hep-ph/0303163 and references therein 

64. Mattingly AC, Stevenson PM. Phys. Rev. D 49:437 (1994) 

65. Dokshitzer YL. In Proc. 29th Int. Conf. on High-Energy Physics (ICHEP 
98), Vancouver, Canada, Jul 23-29, 1998. High energy physics. Vol. 1, p. 305, 
|hep-ph/9812252| 

66. Brodsky SJ. New Perspectives for QCD. . . . |hep-ph/0310289| (2003) 

67. Grunberg G. Phys. Rev. D 29:2315 (1984) 

68. Dokshitzer YL. Light and Heavy Quark Jets in Perturbative QCD. In Physics 
up to 200 TeV, p. 17. Erice (1990); Dokshitzer YL, Khoze VA, Troyan SI. J. 
Phys. G: Nucl. Part. Phys. 17:1481 (1991); Phys. Rev. D 53:89 (1996) 

69. Dokshitzer YL. Why It Seems Too Early to Report on Status of QCD. In 
From Super symmetry to the Origin of Space- Time, p. 108. Erice (1993) 

70. Pomeranchuk IY, Smorodinsky YA. Journ. Fiz. USSR 9:97 (1945); Pieper 
W, Greiner W. Z. Phys. 218:327 (1969); Zeldovich YB, Popov VS. Uspekhi 
Fiz. Nauk 105:4 (1971) 

71. Gribov VN. Possible solution of the problem of quark confinement. Presented 
at Perturbative QCD Workshop, Lund, Sweden, May 21-24, 1991. LU-TP- 
91-7 (unpublished). In ft®, pp. 357-412 

72. Gribov VN. QCD at Large and Short Distances. Eur. Phys. J. C10:70 (1999), 
hep-ph/9807224| 



73. Vacuum Structure and QCD Sum Rules: Reprints, ed. M.A. Shifman. Vol. 



60 Dokshitzer <fe Kharzeev 

10 of Current Physics, Sources and Comments. North-Holland (1992) 



74. Ewerz C. Eur. Phys. J. C13:503 (2000), hep-ph/0001038 



75. Gribov VN. The Theory of Quark Confinement. Eur. Phys. J. C10:91 (1999), 
hep-ph/9902279 



76. Binosi D, Papavassiliou J. J. Phys. G30:203 (2004), hep-ph/0301096; Wat 



son NJ. The Pinch Technique Approach to Gauge Independent N- Point 
Functions. In Proc. Int. Workshop on Physical Variables in Gauge Theories, 



Dubna, Russia, Sep 21-25, 1999, hep-ph/9912303 



77. Gribov VN. Quantum Electrodynamics at Short Distances. (1996, unpub- 
lished). In pp. 519-554 

78. Gribov VN. Higgs and Top Quark Masses in the Standard Model without 
Elementary Higgs Boson. Phys. Lett. B336:243 (1994). In (79), pp. 473-481; 
Bound States of Massless Fermions as a Source for New Physics. In Proc. 
33 rd Int. School of Subnuclear Physics, Erice (1995). In dJ), pp. 483-496. 

79. Gribov VN. Gauge Theories and Quark Confinement. Collection of works. 
Moscow: Phasis Publishing house (2002), 640 pp. 

80. Gribov VN. Orsay lectures on confinement. (I) The Theory of Supercharged 
Nucleus. LPTHE-ORSAY 92-60 (1992). In GS), pp. 421-432; (II) The Con- 
finement of the Heavy Quark. LPTHE-ORSAY 94-20 (1994). In pp. 
433-440; (III) Light Quark Confinement. LPT-ORSAY 99-37 (1999). In (jT^ . 
pp. 441-454 



